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Abstract. - We present a general stochastic dynamic model of a cluster of biological com- 
plexes with fluctuating dissociation and association rates. The master equation has analytical 
solutions in two limiting cases: the vanishing force with reflecting boundary condition and slow 
reaction limit. In particular, the latter can describe the dynamics of a catch-slip bond cluster 
under a constant force. We find that if the rebinding rate vanishes, the size of the cluster mat- 
ters little to the cluster lifetime within the whole catch force regime. The large cluster rapidly 
decays and single catch bond governs the final dynamic behaviors of the cluster. The catch 
behavior would be further amplified in the presence of large nonvanishing rebinding rates. 



Introduction. - Recent experiments found that at single molecule level [1-3] , the dissoci- 
ation rates of some biological adhesion complex, including the bonds between L- and P-selectin 
and P-selectin glycoprotein ligand 1 (PSGL-1) firstly decrease with initial application of small 
force, and then increase when the force increases beyond some threshold. This counterintu- 
itive behavior has been termed catch-slip bond transitions. Because these bonds are primarily 
responsible for the tethering and rolling of leukocytes on inflamed endothelium under shear 
stress [4,5], this discovery may provide a direct experimental evidence at the single-molecule 
level to account for the shear threshold effect [6, 7], in which the number of rolling leukocytes 
first increases and then decreases while monotonically increasing shear stress. Several kinetic 
models [2, 8, 9] have been developed to explain the catch-slip bond behavior at the single 
molecule level. We [10-12] recently suggested alternative mechanism. We proposed that the 
forced dissociations of PSGL-1— L- and — P might a typical rate process with dynamic disor- 
der [13]. The catch-slip bond transitions can be induced by a particular energy barrier shape 
with respect to the inner conformational coordinate [11]. Although applied force projection on 
the dissociation reaction only accelerates the dissociations, the component of the force along 
the conformational coordinate could stabilize the molecular complex by dragging the molecule 
to higher energy barriers; if effect of the latter is greater than that of the former, the bond 
presents catch behavior, otherwise a slip bond is observed 

But before we apply the concept of the catch-slip bond transitions to the real problems, 
we must know that, in physiological setting adhesion complexes usually exist in a cluster and 
cooperatively realize their functions with clusters [14]. Two new features will emerge in a 
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cluster of bonds. Firstly the effect of bond association has to be considered. For single bond 
case, in general this effect could be neglected due to elastic recoil of the force transducer after 
complex rupture [15]. In contrast, the dissociated complex can be formed again when the 
other bonds are still closed. On the other hand, because the closed complexes share the same 
external force, any dissociation or association of one bond would alter the partition of the 
force on the other closed bonds. It means that we must consider the cooperativity between 
bonds in clusters. There are considerable theoretical works contributed to the behaviors of a 
cluster of slip bonds [14, 16-24]. We naturally ask what difference will be observed between 
a cluster of catch-slip bonds and slip bonds. We try to give an answer in the present work. 
Considering that the dissociation of a catch-slip bond is a rate process with dynamic disorder, 
we begin with a general stochastic dynamical model of a cluster of molecular complex with 
fluctuating dissociation and association rates. 

Model. - Following the previous theoretical work [11], we assume that the forced disso- 
ciation rate of single bond has the Bell expression kj{x) = k~^{x) exp(/3/^t), where x is the 
complex conformational coordinate [25], P = UbT, fee is the Boltzmann's constant, T is abso- 
lute temperature, is a projection of the distance from the bound state to the energy barrier 
onto the external force /, k'^{x) is the intrinsic dissociation rate and has Arrhenius form, i.e., 
fc+(x) cx exp [/3AG-'-(x)] . The rebinding rate k~ {x) is independent of force [14, 17,24]. Given 
a cluster of N parallel catch-slip bonds, if the closed bonds sharing the constant force equally, 
the probability that i bonds are closed and has particular conformational coordinate Xi at 
time t then satisfies the following one-step master equation 
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where Ui and bi respectively denote the ruptured and closed bond number, the superscript 
(Tj = 1, respectively represents z-bond ruptured or not, and at = m, 2 < m < N — 1; 
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the Fokker- Planck operators are defined by 
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where D is the difTusion coefRcient, and the potential Vf{x) — Ui{x) — fx [11]. Because 
rebinding of the completely dissociated state is usually prevented by elastic recoil of the 
transducer, we write the first two equations (m=0, 1) independently. According to previous 
theoretical work [24] , they were termed as the absorbing boundary conditions of the master 
equation at m = and 1. If at initial time, all bonds arc closed and in thermal equilibrium, 
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and others vanish, where the distribution P^^{xi) oc exp [—f3Ui{x)], we'd like to calculate the 
cluster lifetime, T/ = tDf{t)dt. Here 
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is an extended cluster dissociation rate with dynamic disorder. We see this master equation are 
far complicated than those in the absence of dynamic disorder. It is impossible to analytically 
solve eq. for general rate functions and potentials with respect to the conformational 
coordinate. Numerical approach such as Monte Carlo simulations [24] is also not simple to 
be applied in the current case. In the following sections, we consider only two limiting cases, 
the vanishing force with reflecting boundary conditions and slow reaction limits, in which the 
master equation has analytical solutions. 

Vanishing force with reflecting boundary conditions. - In this case, the reflecting bound- 
ary condition [24] means that the first two equations in eq. Q are altered into, 
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The other terms (m > 2)are the same with those in eq. (QJ except the force therein vanishing. 
It is easy to prove that Pm{x'[^ ■ • ■ x'^ , t) = -f"^' i^i, t) is the solutio in this limit, where 
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P"'{x,t), a — 0,1 are solutions of the dissociation-association diffusion equations of single 
bond in the absence of force given by 

dP^x,t) 



dt 

dP°{x,t) 



= C} (x)P^ (x, t) - k+ {x)P^ (x, t) + k- (x, t), 

= /:°(x)P°(x,t) + fc+(.T)Fi(x,t)-fc-(.T)F°(x,0. (5) 

Here the initial conditions are P^(a;,0) oc exp(— /3J7i(a:;)) and P'^(.x, 0) = 0. Therefore, any 
quantity of interest can now be calculated. For instance, the probability Qm{t) that m-bond 
are closed at time t is 

N 

Qrnit) = J PmK' ■■■X%",t)l[dx, 

(Tl-'-CTiV 1 

= Cn[ P\x,t)dx] ( / P%x,t)dx] . (6) 

Although the master equation with reflecting boundary conditions has been formally solved, 
it is still a great challenge to solve eq. lO analytically. Various approximation approaches 
have been developed in the past two decades [26]. We are not ready to further discuss it. 
Instead, we study the other limit, in which is relevant to the dynamics of a catch-slip bond 
cluster. 

Slow reaction limit. - In this limit the processes of the rupture and rebinding of any 
single bonds are very slow compared to the conformational diffusions. The thermal equilibrium 
distribution of the conformation for each bond hence is always maintained during the courses 
of "reactions". We assume 

Pm{x^,' ■■■X%-,t)=qra{ai---CTN,t)l[ P^\ (xfc. , //m) J] ^oq i^u^ ) , (7) 
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where P^^, tr = 0, 1 are the thermal equilibrium distributions and satisfy 

4/™(^^.)^eq {=Cb,J/m) = 0, and £°(x„jF°q(a;„J = 0, (8) 

respectively. Substituting in eq. Q and integrating all conformational freedoms lead to a 
master equation about qm{<^i ■ ■ ■ (^n ,t)- According to the definition of Qm{t) (eq. ^),we can 
further simplify the new master equation as follows, 
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We see eq. is formally the same with the one-step master equation that describes the force 
ruptured dynamics of a cluster of slip bond [24] . Many analytical results that are independent 
of concrete expression of the rates hence could be applied here. On the other hand, because 
the new defined rate constant K^{f) is a dynamical average of the Bell expression on the 
distribution Pcq{x, f /m), its force dependence could be very different from that of conventional 
slip bond. Particularly, if the potential Ui{x) is harmonic with spring constant k, and barrier 
AG^{x) is a piecewise function with two segments, K^{f) is approximated to be 



■mkcexp{-l3dlf/m) f < ft 
inks e:xp{+f3dlf/m) f > ft, 



where /( is the bond transition force, ki, i — d,s are some constants, and the "effective" 
distance dj > 0, are the simple combinations of k, the slopes of segments, and [11]. We 
see if force is smaller than ft, the rupture rate K^^{f) decreases with increasing force (a catch 
bond cluster); otherwise it is a shp bond cluster. For the latter, eq. © is completely the same 
with previous models [22,24]. Hence we only study the dynamics of the catch bond cluster. 
We respectively discuss two cases according to zero and nonzero rebinding rates. 

Rebinding rates k~{x) = 0. In this case, the probability Qm{t) of m-bond closed at time 
t has an exact solution given by [22] 
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Hence, any quantity of interest can be calculated from Qm{t)- For instance, the cluster lifetime 
Tf is 
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For large N and small force f < fc — kBT/dc, the above equation can be approximated to 
be [28] 

T/ « i/jv « In N + 1/2N + 0.577, (13) 

where Hjy — harmonic number. It is same with the lifetime of slip bond 

cluster [24]. It is expected because at small forces the difference between the catch and slip 
bonds are negligible. If / > fc, the force in the exponentials of Tf cannot be neglected again. 
In contrast, for the larger force the dominating term in eq. H12f) is the one for the rupture of 
only one bond, i.e., 

Tf^kc^'eMf/fc). (14) 

This result is further demonstrated by the time decay function of the mean number of the 
closed bonds on large time scales: 

(m) ~ exp [-tkc exp{-f/fc)] ■ (15) 



Compared to the behaviors of a slip bond cluster with large size, three difference are observed 
in a large size of catch bond cluster if the rebinding rate vanishes. First, the size of the 
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catch bond cluster matters little for the cluster lifetime within the whole catch force regime, 
whereas it only holds for smaller forces in the slip bond cluster. Then, for larger force the 
dynamics of the catch cluster is took over by single catch bond. In contrast, the largest clusters 
dominate the decay of the slip bond cluster [24] . Finally, as a characteristic of a catch bond, 
the increasing force only exponentially elongates the lifetime of the catch bond cluster. 

Nonzero rebinding rates. Although Qm{t) in Eq. ^ in principle could be calculated by 
eigenvalue technique in this case, except the simplest case A'^ = 2, it is not more efficient than 
the exact stochastic simulations [24]. Fortunately, it has been proved [24,27] that the cluster 
lifetime of great interest has still an exact form 



We are not interested in the general dependence of eq. H16() on force and N. Instead we 
consider a specific case where the ratio 7 = {k^)/kc is larger than 1. Thus for large N and 
force, eq. lfTC|) is 



In contrast to the case without rebinding, the size of the catch cluster will change the cluster 
lifetime dramatically as 7 > 1; the presence of the rebinding rate further amplifies the lifetime 
by 7^~-'^-fold. Therefore, the effect of the rebinding are the same for the slip and catch bond 
cluster. 

Discussion and conclusions. - In this work, we present a general stochastic dynamic 
model of a cluster of biological complexes with dynamic disorder. We study two simplest 
limiting cases: vanishing force with reflecting boundary conditions and the slow reaction 
limits. Under the former limit, because the dissociation and association of each bonds are 
independent, the solution to the model is a simple linear combination of the dynamics of 
single bond. Under the slow reaction limit, we demonstrate that the model can be reduced 
to the conventional one-step master equation except that the rates are replaced by dynamic 
averages. Particularly, given a specific energy barrier with respect to molecular conformational 
coordinate, this master equation can describe the dynamics of a cluster of catch-slip bonds. 
We found that many conclusions obtained for the slip bond cluster could be applied to the 
catch bond cluster with minor modification. We give simple relations between the cluster 
lifetime and the cluster size N, force / and the average rebinding rate. These conclusion 
would be useful in future dynamic studies of the catch-slip bond clusters. 
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